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Abstract 

■ We deal with the Hamiltonian hierarchy problem of the Hulthen poten- 

tial within the frame of the super symmetric quantum mechanics and find 
| that the associated superymmetric partner potentials simulate the effect 

\Q ' of the centrifugal barrier. Incorporating the super symmetric solutions and 

using the first-order perturbation theory we obtain an expression for the 
energy levels of theHulthen potential which gives satisfactory values for 
the non-zero angular momentum states. 
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1 Introduction 

TheHulthen potential [1,2] is one of the important short-range potentials in 
physics. The potential is given by 

where Z is a constant and S is the screening parameter. If the potential is used 
for atoms, the Z is identified with the atomic number. This potential is a special 
case of the Eckart potential || which has been widely used in several branches 
of physics and its bound-state and scattering properties have been investigated 
by a variety of techniques (see e.g., j|] and references therein). 

Unfortunately, the radial Schrodinger equation for the Hulthen potential can 
be solved analytically only for the states with zero angular momentum [1,2,5,6]. 
For i ^ 0, a number of methods have been employed to evaluate bound-state 
energies numerically [4,7-17]. In connection with this, we present in this letter 
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a method within the frame of supersymmetric quantum mechanics (SUSYQM) 
using an effective Hulthen potential for non-zero angular momentum states, which 
can be solved analytically. 

In supersymmetric quantum mechanics (for a recent review see ||18|| ) one of- 
ten deals with hierarchy problem. Within the context of the SUSYQM one can 
generate a Hamiltonian hierarchy, the adjacent members of which are supersym- 
metric partners in that they share the same eigenvalue spectrum except for the 
missing ground state. In the case of Coulomb potential V c (r) , the Hamiltonian 
hierarchy corresponds to the addition of an appropriate centrifugal potential and 
the so-called accidental degenarcy is recovered as a natural consequence |19|]. In 



this letter we shall examine the implication of the Hamiltonian hierarchy for the 
Hulthen potential. At small values of the radial coordinate r, the Hulthen poten- 
tial behaves like a Coulomb potential whereas for large values of r it decreases 
exponentially so that its capacity for bound state is smaller then V c (r). In con- 
trast to the Hulthen potential, the Coulomb problem is analytically solvable for 
all energies and all angular momenta. Because of the similarity and points of 
contrast mentioned above between Coulomb and Hulthen potentials, it may be 
of considerable interest to generate the supersymmetric partners of the Hulthn 
potential and study their eigenf unctions and eigenvalues. In the following we 
outline the basic idea of the SUSYQM and set up the formalism for dealing with 
the Hulthen problem. 



2 Supersymmetric solution of theHulthen po- 
tential 

The supersymmetric formalism has already been used to study some aspects of 
the Hulthen potential [15,17]. Here the exact analytical solution for this potential 
is re-obtained for £ = state in the light of the works described in Refs. [15- 
17,20] to show the consistency of the method and see how partners of the Hulthen 
potential simulate the effect of the centrifugal barrier, which makes clear the 
physics behind partner Hamiltonians in connection with the states having £ ^ 
angular momenta. 

In the approach followed here the first step taken is to look for an effective 
potential similar to the original Hulthen potential. Inspired by SUSYQM, we 
propose an ansatz for the superpotential, 



m H (£+l)6e~ Sr [me 2 

W t+ i(r) = —^=-7-, TTTT" + J TT 

V2m (1-e drj V 2 h 



1 (£ + l)P 
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where (^+1) denotes the partner number with £ = 0, 1, 2, and (3 = ^| which is 
a dimensionless quantity. This kind of superpotential choice leads to the (£+l)-th 
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member of the Hamiltonian hierarchy: 
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We introduce an expression for the bound-state energies of the above potential, 



considering the shape invariance requirement |2T 



E, 



me 



2h z 
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and the corresponding ground-state eigenfunctions are 
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It is reminded that for a number of purposes it is convenient to have the wave- 
function in such a compact analytical form. The first eigenfunction corresponds 
to the minimum energy for each I. In terms of the hierarchy of Hamiltonians, we 
present here the lowest state wavefunctions for each member. The excited state 
wavefunctions can be determined [18,26] from the usual approach in SUSYQM. 

For t = the potential in Eq. (4) leads to the usual Hulthen potential which 
has an interesting property such that when the angular momentum is zero it is 
not shape invariant in the sense expressed in |^TJ. However, it is still possible 
to construct a general form of the potentials in the shape invariant super-family 
of Hamiltonians as seen in Eq. (4) where the first member corresponds to the 
Hulthen potential which can be solved exactly in analytic form. One can easily 
verify that the energy eigenvalues and eigenfunctions for I = case of Eqs. 
(5) and (6) are the same given in Refs. [5,6]. This supports the suggestion 
[17,22-24] that the Gedenshtein's condition of shape invariance is sufficient but 
not a necessary condition in the construction of exactly solvable but non-shape 
invariant potentials. 

Eq.(4) can be rearranged as 



Vg f (r) = V i+1 (r) 



5e 
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1-e 
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which is known in literature as the approximate Hulthen effective potential in- 
troduced by Greene and Aldrich El in their method to generate pseudo-Hulthen 
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wave functions for £ ^ states. For small Sr, Eq. (7) is a good approximation 
to the realisticHulthen effective potential, and unlike the original case the radial 
Schrodinger equation for this potential is solvable analytically through Eqs. (5) 
and (6). In addition, the partner potentials in Eq. (7) gives the necessary repul- 
sive core due to angular momentum. For instance, for small r the second term 
in Eq. (7) behaves as a p-wave centrifugal barrier for the second member of the 
super-family. Since we know that the centrifigal potential is effective only in this 
region (i.e., small r), eigensolution of the potential for i — 1 in Eq. (6) can be 
regarded as the approximate £>-wave solution for the Hulthen potential. Clearly, 
one can get other supersymmetric partners and their solutions in an explicit form 
for £ ^ states. The present simple and elegant method is a clear cut of the 
iteration technique introduced by Laha et al. [15,16] 

For the sake of completeness, it is of interest to note that for small values of 
S, the potential in Eq. (7) closely approximates the effective Coulomb potential 
rather well, 



Vg f {r, 5 = 0)^ Vc ff {r) = -— + 



£{£+l)h 2 
2mr 2 



(8) 



and the corresponding energy eigenvalue for the potential of Eq. (8), together 
with its ground state wavefunction for £ ^ states, obtained easily via Eqs. 
(5) and (6) overlap with those, e.g. in Ref. [ 26| . This makes clear the work 
of Lam and Varshni H in which they showed that if one uses as trial functions 
eigenvectors of theHulthen potential rather than those of the simple Coulomb 
potential, excellent results for the energies of the states of the screened Coulomb 
potential can be obtained with simple variational wave functions containing only 
one parameter. 

An important quantity of interest for the Hulthen potential (and for other 
similar screened potentials) is the critical screening parameter 5 C , which is that 
value of 5 for which the binding energy of the level in question becomes zero. 
Using Eq. (5), in atomic units, 



(n + e + iy 



(9) 



which works well for all n values in case £ = when compared to those in Table 
III of Ref. [|J], but fails for non-zero angular momentum states. Consequently, 
the eigenenergies obtained via Eq. (5) for £ ^ states deviates from the accurate 
values obtained by numerical techniques and presented in Table I of Ref. H. 
This may be understood as follow. If Eq. (7) is written in the form 
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the exact energy eigenvalues for the realistic effective Hulthen potential may be 
given as 



E ni _ me 



1 {n + £+l)(3 l2 



(n + 



+ AE (11) 



where AE is the contribution, which does not appear in Eq. (9), due to the 
last term in Eq. (10). The clear interpretation of Eqs. (10) and (11) is that 
the potential barrier term prevents us to build the super-family as in the £ = 
case, since the potential-the first two terms in Eq. (10)-is not exactly solvable 
hence the supersymmetry is broken for £ ^ due to the potential barrier term. 
It is easy however to verify that for small values of 5, AE goes to zero while Eq. 
(10) becomes an expression for the effective Coulomb potential in which case the 
accidental degenarcy is recovered as a natural consequence. 

The usefulness of theHulthen potential would be enhanced if one obtains an 
analytical expression for the exact energies of the non-zero angular momentum 
states. The work along this line is in progress in the frame of broken supersymetry. 
Further, in the light of the supersymmetric solutions discussed in this letter we 
suggest here, as an alternative to other various methods [4,7-14] investigating 
the bound-state properties of theHulthen potential, an elegant approach for the 
calculation of the whole energy spectrum of the potential using the first-order 
perturbation theory, 



E H = E e+1 + 



2m 



I" 00 r 1 2 ( 1 ffi \ 



which gives satisfactory values when compared (see Table 1) with the results 
obtained by the various methods for the eigenenergies of £ ^ levels. The 
accuracy of the present calculations may be improved incorparating higher-order 
perturbations for in particular large values of the screening parameter. 



3 Conclusion 

We have obtained the exact analytical eigenfuntions and eigenvalues for the 
Hulthen potential within the framework of SUSYQM for the case £ = 0. The 
approach consists of making an ansatz in the superpotential which satisfies the 
Riccati equation by an effective potential. For £ = the effective potential ob- 
tained is identical to the Hulthen potential. However, for £ ^ the effective 
supersymmetric potential has a slightly different structure than the Hulthen po- 
tential. This deviation has led us to introduce a simple expression that yields 
reasonable results for the non-zero angular momentum state energies. We stress 
that even though the problem has been attacked by different methods our simple 
and elegant methodology is powerful because it provides an insight into the re- 
lation between theoretical partner Hamiltonians in the frame of supersymmetric 
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quantum mechanics and physical states of the system considered. We hope to 
stimulate further examples of applications of SUSYQM in important problems of 
nuclear and atomic physics. 
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Table 1: Energy eigenvalues of the Hulthen potential as a function of screening 
parameter for various states in atomic units. 



State 5 Present Numerical Variational Lai Patil Tang Matthys 

Calculations Integration (Ref. [4]) and Lin (Ref. [9]) and Chan and De Meyei 
(Ref. [4]) (Ref. [7]) (Ref. [12]) (Ref. [14]) 



2p 


0. 


.025 


0.1127605 


0.1127605 





.1127605 






0.11276 







.1127604 




0. 


.050 


0.1010425 


0.1010425 


0. 


.1010425 





.101043 


0.10104 


0.1010424 





.1010425 




0. 


.075 


0.0898478 


0.0898478 


0. 


.0898478 






0.08985 










0. 


100 


0.0791794 


0.0791794 


0. 


.0791794 





.079179 


0.07918 


0.0791794 





.0791794 




0. 


.150 


0.0594415 


0.0594415 


0. 


.0594415 






0.059445 







.0594415 




0. 


.200 


0.0418854 


0.0418860 





.0418860 





.041886 


0.041895 


0. 0418857 





.0418860 




0. 


.250 


0.0266060 


0.0266111 


0. 


.0266108 
















0. 


300 


0.0137596 


0.0137900 


0. 


,0137878 


0. 


.013790 









.0137900 




0. 


350 


0.0036146 


0.0037931 





.0037734 





.003779 


0.038375 








3p 


0. 


.025 


0.0437068 


0.0437069 


0. 


.0437069 





.043707 


0.0437085 







.0437071 




0. 


050 


0.0331632 


/ \ / \ i k k 1 1 / * At" 

0.0331645 


0. 


.0331645 





.033165 


0.033185 


0.03316518 


0. 


.0331650 




0. 


.075 


0.0239331 


0.0239397 





.0239397 






0.0240165 










0. 


100 


0.0160326 


0.0160537 





.0160537 





.016054 


0.01622 


0.01606772 


0. 


.0160537 




0. 


150 


0.0043599 


0.0044663 





.0044660 





,004466 


0.046995 




0. 


.0044664 


3d 


0. 


.025 


0.0436030 


0.0436030 


0. 


,0436030 





.043603 


0.0436025 







.0436030 




0. 


.050 


0.0327532 


0.0327532 





.0327532 





.032753 


0.032745 


0.0327532 





.0327532 




0. 


.075 


0.0230306 


0.0230307 





.0230307 






0.02299 










0. 


100 


0.00144832 


0.0144842 


0. 


.0144832 





.014484 


0.01439 


0.0144842 





.0144842 




0. 


150 


0.0132820 


0.0013966 


0. 


.0013894 





.001391 


0.0013755 







.0013965 


4 P 


0. 


.025 


0.0199480 


0.0199489 





.0199489 





.019949 


0.01995 







.0199490 




0. 


050 


0.0110430 


0.0110582 


0. 


,0110582 


0. 


.011058 


0.011075 


0.0110725 





.0110583 




0. 


.075 


0.0045385 


0.0046219 


0. 


,0046219 


0. 


.004622 


0.0046585 







.0046224 




0. 


100 


0.0004434 


0.0007550 


0. 


,0007532 


0. 


.000754 


0.000752 








4d 


0. 


.025 


0.0198460 


0.0198462 





.0198462 





.019846 


0.019845 







.0198462 




0. 


050 


0.0106609 


0.0106674 


0. 


.0106674 





.010667 


0.01068 


0.0106690 





.0106674 




0. 


.075 


0.0037916 


0.0038345 


0. 


.0038344 





.003834 


0.003875 







.0038346 


4f 


0. 


.025 


0.0196911 


0.0196911 





.0196911 





.019691 


0.01969 







.0196911 




0. 


050 


0.0100618 


0.0100620 


0. 


,0100620 





.010062 


0.010045 


0.0100620 





.0100619 




0. 


.075 


0.0025468 


0.0025563 


0. 


,0025557 





.002556 


0.002557 







.0025563 
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Table 2: Continue of Table 1. 



State 5 Present Numerical Variational Lai Patil Tang Matthys 

Calculations Integration (Ref. [4]) and Lin (Ref. [9]) and Chan and De Meyer 
(Ref. [4]) (Ref. [7]) (Ref. [12]) (Ref. [14]) 



5p 0.025 0.0094011 0.0094036 
0.050 0.0026056 0.0026490 



0.0094087 



5d 0.025 0.0092977 0.0093037 
0.050 0.0022044 0.0023131 



0.0093050 



5f 0.025 0.0091507 0.0091521 
0.050 0.0017421 0.0017835 



0.0091523 



5g 0.025 0.0089465 0.0089465 
0.050 0.0010664 0.0010159 



0.0089465 



6p 0.025 0.0041493 0.0041548 

6d 0.025 0.0040452 0.0040606 

6f 0.025 0.0038901 0.0039168 

6g 0.025 0.0036943 0.0037201 
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